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►  A  new  method  for  creating  automated  parameterized  time-domain  models  is  introduced. 

►  Impedance  measurements  in  a  broad  frequency  range  are  performed. 

►  The  Distribution  of  Relaxation  Times  (DRT)  is  calculated. 

►  The  model  is  successfully  validated  with  different  current  profiles  in  the  time-domain. 

►  The  advantages  and  limits  of  this  modeling  approach  are  discussed. 
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A  novel  approach  for  the  development  of  a  model  for  Li-ion  batteries  with  the  potential  for  application  in 
on-board  diagnostic  is  introduced.  This  approach  is  not  tied  to  a  particular  electrochemistry,  and  has  the 
advantages  of  scalability  and  automatable  parameterization  as  well  as  feasibility  regarding  an  imple¬ 
mentation  on  microcontrollers.  The  dynamic  behavior  of  the  cell  under  test,  a  Li-ion  pouch  cell  with 
a  graphite  anode  and  a  blend  cathode,  is  characterized  in  an  extremely  broad  frequency  range  from 
100  kHz  to  17  pHz.  The  introduced  model  gives  predictions  consistent  with  the  experimental  data,  its 
accuracy  depends  on  the  assumed  number  of  relaxation  times  respectively  model  order.  Recommen¬ 
dations  concerning  application  are  given  and  limitations  of  the  approach  are  discussed. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

There  are  special  requirements  for  automotive  application  of 
battery  systems.  This  has  led  to  an  increased  number  of  publications 
covering  methods  for  simulation  and  modeling  of  lithium  ion 
batteries.  The  range  of  complexity  and  parameterization  effort  is  very 
broad.  Our  approach  features  the  distribution  of  relaxation  times  and 
is  thereby  bridging  the  gap  between  physically  based  mathematical 
models  and  coarse  equivalent  circuit  models.  It  combines  straight¬ 
forward  parameterization  and  a  structure  that  is  readily  imple- 
mentable  on  microcontrollers  with  a  scalable  degree  of  detail. 

The  application  of  models  is  always  associated  with  constraints, 
i.e.  accuracy,  computational  complexity,  configuration  effort  or 
analytical  insight  [  1  ] .  As  there  is  no  “correct”  model,  the  goal  must  be 
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to  select  the  model,  which  is  optimal  under  given  constraints.  In 
system  theory  models  are  ranked  in  order  of  their  degree  of  physical 
interpretability:  white  box,  gray  box  and  black  box  models  [2]. 

White  box  models  are  based  on  physical  and  electrochemical 
partial  differential  equations.  They  describe  transport  processes 
and  electrochemical  processes,  which  can  be  linked.  These  models 
are  often  simulated  applying  FEM-Tools,  whereby  spatial  resolution 
is  variable.  Calculation  effort  increases  from  OD-models  over 
homogenized  or  particle  ID  models  [3,4]  to  homogenized  2D 
models  [5,6]  up  to  3D  models  [7].  Material  constants  for  parame¬ 
terization  are  either  available  from  literature  or  obtained  from 
extensive  measurements. 

Equivalent  circuit  models,  for  example,  belong  to  the  class  of 
gray  box  models.  They  are  proposed  in  variable  degrees  of  detail. 
Physical  insight  can  be  gained  by  using  adequate  circuit  elements, 
i.e.  for  electrical  conductivity,  charge  transfer,  electrochemical 
double  layer  or  diffusion  [8].  However,  ambiguity  of  equivalent 
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circuit  models  is  a  problem.  Despite  correct  description  of  the 
clamp  behavior,  physical  interpretation  of  circuit  element  values 
can  be  misleading  [9].  In  some  cases  circuit  elements  are  of  frac¬ 
tional  nature  and  cannot  be  simulated  directly  in  time  domain.  This 
can  be  solved  by  the  implementation  of  a  fractional  integrator  [10- 
12]  or  by  approximating  the  system  with  a  sum  of  RC  elements 
[8,13].  While  generally  the  spatial  resolution  of  electrode  processes 
in  equivalent  circuit  models  is  lower  compared  to  white  box 
models,  still  ladder  type  models  can  be  used  to  describe  porosity 
[14-17].  A  coarse  spatial  resolution  is  achieved  via  segmentation,  as 
proposed  in  Refs.  [18,19].  In  general,  equivalent  circuit  models  are 
parameterized  by  fitting  them  to  a  measured  electrochemical 
impedance  spectrum  [13,17,20-25].  Alternatively,  a  current  or 
potential  profile  is  used  for  parameterization  [26-30]. 

Artificial  neuronal  networks  [31—34],  fuzzy  logic  [35]  or  NAR- 
MAX  models  [36]  are  rated  among  black  box  models ,  which  are 
mostly  used  for  diagnosis.  Certainly  their  parameters  do  not  allow 
for  a  physical  interpretation,  but  on  the  other  hand  parameteriza¬ 
tion  can  be  completely  automated  and  no  a  priori  knowledge  of  the 
processes  in  the  batteries  is  needed. 

The  introduced  model  combines  the  advantages  concerning 
physical  insight  of  a  gray  box  model  with  the  fully  automatable 
parameterization  of  a  black  box  model.  It  is  based  on  the  idea,  that 
the  impedance  of  a  capacitive  electrochemical  system  can  be 
interpreted  as  a  continuous  distribution  of  RC  elements  in  the  space 
of  relaxation  times.  This  is  a  very  general  assumption,  which  is  not 
tied  to  a  particular  electrochemistry:  Batteries  and  even  fuel  cells 
have  proven  to  be  well  described  by  models  based  on  that  approach 
[37-40].  Its  distinctive  property,  the  ability  to  nicely  visualize 
physical  properties  of  the  studied  system,  is  based  on  the  improved 
separation  of  physical  processes  like  charge  transfer  or  diffusion  in 
the  space  of  relaxation  times  [41]. 

There  every  process  is  represented  as  local  maximum  in 
a  continuous  distribution  function.  However,  one  does  usually  not 
operate  with  a  continuous  distribution  function.  Rather,  a  dis¬ 
cretization  of  this  function  is  used,  assuming  the  impedance  to 
consist  not  of  an  infinite,  but  only  a  finite  number  of  RC  elements. 

Choosing  this  number  rather  small,  the  physical  interpretability 
of  the  continuous  distribution  function  can  only  be  preserved  by 
applying  a  priori  knowledge  of  the  time  constants  of  the  physical 
processes,  compare  Fig.  lb)  and  c).  Without  a  priori  knowledge,  the 
physical  interpretability  of  the  model  can  only  be  preserved  by 
using  a  sufficiently  large  number  of  evenly  distributed  RC  elements. 
For  the  sake  of  a  fully  automated  parameterization  process,  in  this 
paper  the  latter  method  will  be  described. 

In  the  following  sections  (i)  scalability  and  (ii)  automatable 
parameterization  will  be  discussed.  After  the  experimental  chapter 
the  model  is  validated  with  different  current  profiles  and  limita¬ 
tions  of  the  model  are  discussed. 


2.  Theory 

The  method  of  the  distribution  of  relaxation  times  has  been 
successfully  applied  for  years  for  the  analysis  of  impedance  spectra 
of  solid  oxide  fuel  cells  (SOFC)  [41,42]  and  has  recently  been  applied 
for  the  analysis  of  lithium  ion  cells  [43,44].  Its  major  benefit  is  the 
better  separation  of  processes  with  different  time  constants 
compared  to  the  Nyquist  or  Bode  plot  [41  ].  An  impedance  spectrum 
can  be  written  as  integral  equation,  where  g(i)  represents  the 
distribution  of  relaxation  times: 

oo 

Z(W)  =  R0  +  RpoXj  ^-Ax.  (1) 

0 


f  [Hz] 


Fig.  1.  a)  Continuous  distribution  of  relaxation  time  for  an  exemplary  system  exhib¬ 
iting  three  different  time  constants  and  an  approximation  with  a  numerical  calculated 
discrete  distribution  with  an  evenly  resolution  of  50  time  constants,  b)  The  continuous 
distribution  in  comparison  with  a  discrete  distribution  comprising  only  four  evenly 
distributed  time  constants,  c)  The  continuous  distribution  function  approximated  with 
three  RC-elements  whose  time  constants  were  selected  using  a  priori  knowledge. 


Since  the  inversion  is  performed  numerically,  the  distribution  itself 
is  not  obtained  as  an  analytical  function  but  as  discrete  values. 
Therefore  the  integral  equation  becomes  a  sum: 

Z(w)  =  Ro  +  Rpo,  £  (2) 

Here  <5T>i0g  is  the  difference  between  the  logarithmically  equidistant 
distributed  time  constants  t.  From  this  equation  it  can  be  deduced, 
that  the  measured  impedance  Z(w)  can  be  represented  by  a  serial 
connection  of  N  RC  elements  and  an  ohmic  resistance  Ro-  With  an 
increasing  number  of  relaxation  times  N,  the  number  of  RC 
elements  as  well  as  the  computational  effort  increase. 

The  according  equivalent  circuit  is  depicted  in  Fig.  2.  From 
equation  (2),  the  calculation  of  the  parameters  of  the  nth  RC  element 
can  be  derived.  The  resistance  of  the  RC  element  is  calculated  as 

—  §n  '  ^i,log  ’^pol  (3) 

and  the  capacitance  as 

Gi  —  Tn/^n- 


(4) 
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Fig.  2.  Introduced  equivalent  circuit  model,  based  on  the  distribution  of  relaxation  times  and  extended  for  an  ohmic  resistance  and  a  voltage  source  for  the  description  of  the  open 
circuit  voltage. 


In  this  form  the  model  is  limited  to  capacitive  systems,  whereas 
purely  capacitive  and  inductive  behavior  cannot  be  described. 
Purely  capacitive  behavior  is  observable  for  frequencies 
approaching  zero.  It  is  caused,  for  example,  by  the  differential 
intercalation  capacity  of  battery  electrodes  [45]  and  can  be 
modeled  by  a  voltage  source  (compare  Fig.  2)  [46],  which  repre¬ 
sents  the  open  circuit  voltage  (OCV)  depending  on  SOC. 

If  the  measured  impedance  spectrum  shows  inductive  behavior, 
a  serial  inductance  or  an  RL  element  [47]  has  to  be  fitted  and 
subtracted  prior  to  the  calculation  of  the  DRT.  Afterward  this 
inductance  can  be  added  as  a  serial  element  to  the  equivalent 
circuit  model.  If  the  sampling  rate  of  current  and  voltage  is  below 
100  Hz,  inductive  behavior  cannot  be  observed.  In  this  case  the  last 
step  can  be  omitted. 

For  implementation  the  number  of  relaxation  times  and  the 
covered  frequency  range  have  to  be  chosen  prior  to  the  calculation 


of  the  DRT.  This  allows  for  easy  adjustment  of  the  model  regarding 
complexity,  accuracy  and  calculation  effort.  Since  the  time 
constants  are  logarithmically  equidistantly  distributed,  each 
frequency  decade  is  described  with  the  same  model  order  and  with 
the  same  accuracy.  Hence,  the  dynamic  behavior  is  uniformly 
modeled  for  the  entire  frequency  range  covered.  To  demonstrate 
the  straight-forward  scalability  and  the  resulting  influence  on 
simulation  accuracy  and  speed,  models  of  the  order  N  =  100,  20, 10 
and  4  are  parameterized  and  compared. 

The  procedure  for  parameterization  and  validation  as  well  as  the 
model  structure  is  summarized  in  Fig.  3.  It  contains  the  following 
steps: 

1)  Measurements:  (i)  C/40  discharge  as  quasi-stationary  OCV- 
curve,  (ii)  EIS  measurements  with  a  variation  of  the  SOC,  (iii) 
current  pulses  followed  by  a  relaxation  phase  with  a  variation 


Current 


Scalability 


1^1— 

Ri  R2  Rn  Rn 


Cell  Potential 


Fig.  3.  Flow  chart  showing  the  procedure  of  the  parameterization  of  the  DRT-model.  The  measurement  of  the  validation  profile  is  not  necessary  for  the  parameterization  but  is  used 
for  the  validation  of  the  model. 
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of  the  SOC  for  parameterization  and  (iv)  current  profiles  for 
validation. 

2)  Parameterization:  In  order  to  reduce  measurement  time,  the 
low  frequency  part  of  the  impedance  spectra  is  obtained  by 
transforming  the  current  pulse  measurements  into  frequency 
domain  [48].  Then  impedance  spectra  from  EIS  and  pulse 
measurements  are  merged  and  the  DRT  is  calculated  for  each 
SOC. 

3)  Model:  The  values  Rn  and  Cn  obtained  from  the  DRT  are  used  for 
the  parameterization  of  the  RC-model  implemented  in  Matlab / 
PLECS.  The  quasi-stationary  OCV-curve  from  step  1 )  is  repre¬ 
sented  by  a  controlled  voltage  source. 

4)  Simulation  and  validation:  The  current  profiles  applied  in  step 
1 )  are  now  simulated  and  the  resulting  potential  is  compared 
with  the  measured  potential. 

All  single  steps  are  further  explained  and  discussed  in  the 
section  Results  and  discussion. 

3.  Experimental 

The  investigated  cell  is  a  Li-ion  pouch  cell  from  KOKAM 
(SLPB834374H)  with  a  nominal  capacity  of  2  Ah.  The  anode  consists 
of  graphite;  the  cathode  is  a  blend  of  LiNio.8Coo.15Alo.05O2/LiCoO2. 
All  measurements  were  conducted  in  a  climate  chamber  at  23  °C. 
The  cell  was  contacted  with  a  clamping  device  consisting  of  copper 
blocks,  which  realize  four-terminal  sensing  and  also  serve  as  heat 
sinks.  All  electrochemical  experiments  were  conducted  with 
a  1470E  cell  test  system  (Solartron,  UK).  To  increase  measurement 
accuracy  for  the  OCV-curve,  the  potential  was  additionally 
measured  with  a  34902A  data  acquisition  unit  (Agilent,  UK). 

At  the  beginning,  the  actual  capacity  was  measured  and  the 
stability  of  the  cell  was  verified  by  tenfold  repetition  of  the  capacity 
test.  The  cell  was  charged  with  a  current  of  1C  and  then  the 
potential  was  held  at  4.2  V  for  1  h.  Immediately  hereafter,  the  cell 
was  discharged  with  a  current  of  1C  until  the  stop  voltage  of  2.7  V 
was  reached. 

Next  the  quasi-stationary  OCV-curve  was  recorded:  First  the  cell 
was  charged  with  the  aforementioned  charge  protocol  and  then 
discharged  with  a  current  of  C/40  until  the  stop  voltage  was 
reached.  Then  the  direction  of  the  current  was  reversed  and  the  cell 
was  charged  with  C/40  until  the  upper  stop  voltage  of  4.2  V  was 
reached.  After  integration  of  the  current  and  normalization  using 
the  actual  capacity  from  the  capacity  test,  the  OCV-curve  can  be 
displayed  as  a  function  of  the  SOC,  see  Fig.  4. 

Afterward  impedance  spectra  from  an  SOC  of  95%  down  to  0% 
were  recorded  subsequently  in  5%  decrements.  The  potential 
perturbation  amplitude  was  chosen  at  5  mV  and  the  impedance 
was  recorded  for  frequencies  ranging  from  100  kHz  to  5  mHz  with 
twelve  points  per  frequency  decade.  After  each  EIS  measurement 
a  relaxation  phase  of  5  h  was  inserted,  before  a  current  pulse  of 
-2  A  with  a  duration  of  10  s  was  applied.  The  relaxation  of  the  cell 
potential  was  recorded  for  40  h.  According  to  the  procedure 
described  in  Ref.  [48],  the  measurement  results  from  the  current 
pulse  measurement  have  been  transformed  into  frequency  domain. 
The  applied  algorithm  is  based  on  the  Fourier  transform,  allowing 
the  calculation  to  be  performed  without  prior  selection  of  a  model. 
The  impedance  spectra  from  the  time-domain  measurements  are 
then  merged  with  the  impedance  spectra  from  the  EIS  measure¬ 
ments.  Finally  the  DRT  is  computed  from  the  merged  impedance 
spectra  for  a  frequency  range  of  100  kHz  down  to  as  low  as  17  pHz. 

To  validate  the  model  in  the  time  domain,  after  fully  charging 
the  cell  it  was  discharged  by  repeating  the  current  profile  shown  in 
Fig.  5.  The  maximum  current  Jmax  was  varied  between  C/4  and  4C 
for  different  runs,  which  results  in  an  average  discharge  current  of 


SOC  [%] 

Fig.  4.  OCV  as  function  of  SOC,  measured  at  constant  current  discharge  and  charge 
with  C/40. 

Jmax/2.  In  every  run,  the  discharge  was  stopped  as  soon  as  the  stop 
voltage  of  2.7  V  was  reached.  This  way  the  dynamic  behavior  of  the 
cell  over  the  full  SOC-range  is  evaluated. 

4.  Results  and  discussion 

4.1.  Evaluation  of  model  quality  infrequency  domain 

The  merged  impedance  spectra  from  EIS  and  time-domain 
measurements  are  depicted  in  Fig.  6.  For  a  decreasing  SOC,  the 
polarization  resistance  is  increasing.  The  DRT  can  now  be  calculated 
from  the  obtained  impedance  data.  Prior  to  the  computation  the 
number  of  time  constants  has  to  be  selected.  For  a  number  of 
relaxation  times  N  =  100  and  a  frequency  range  from  3.8  kHz  to 
17  pHz  the  resulting  DRT  is  depicted  in  Fig.  7. 

The  contributions  of  the  polarization  processes  are  distributed 
over  the  whole  frequency  range.  This  overall  distribution  cannot  be 
represented  by  a  standard  gray  box  model  like  for  example  two 
serially  connected  RC  elements.  While  a  dependence  on  the  SOC 
can  hardly  be  stated  for  frequencies  above  10  Hz,  for  frequencies 
between  10  Hz  and  1  Hz  a  distinctive  dependence  is  noticeable. 
Here  the  peak  frequency  decreases  monotonously  with  decreasing 
SOC  and  the  area  beneath  the  peak  increases.  This  behavior  can  be 
attributed  to  the  charge  transfer  processes,  as  it  was  described  in 
a  recent  publication  [49]. 

The  sudden  increase  of  the  DRT  at  low  frequencies  at  an  SOC  of 
0%  correlates  very  well  with  our  observation  that  relaxation  is 
extremely  slow  at  the  end  of  a  complete  discharge.  However,  unlike 


Fig.  5.  Single  cycle  of  the  current  profile  for  the  validation  of  the  model  in  time 
domain.  The  cycles  are  repeated  until  the  cell  potential  reaches  2.7  V. 
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Fig.  6.  Measured  impedance  spectra  at  different  SOCs  with  frequencies  ranging  from 
100  kHz  down  to  17  pHz.  The  frequency  range  from  100  kHz  to  5  mHz  was  investigated 
applying  EIS  measurements,  below  10  mHz  the  impedance  was  investigated  using 
pulse  tests.  The  overlapping  frequency  range  was  used  to  merge  the  spectra  at  10  mHz 
(black  circles),  resulting  in  one  valid  impedance  spectrum. 

the  charge  transfer  process,  in  the  low  frequency  region  a  non- 
monotonous  behavior  is  observed.  Although  this  observation  is 
qualitatively  correct,  it  has  to  be  mentioned  that  the  relaxation 
phase  was  not  long  enough  to  contain  the  highest  relaxation  times, 
especially  for  SOCs  below  30%.  Therefore,  some  peaks  of  the  lowest 
frequency  process  seem  to  occur  at  the  same  lowest  observable 
frequency. 

While  a  physical  interpretation  is  undeniably  interesting  and 
the  DRT  can  be  applied  here  with  the  benefit  of  separation  and 
identification  of  processes,  this  is  not  necessary  for  the  parame¬ 
terization  of  the  model.  The  quality  of  the  parameterized  model  can 
be  assessed  by  comparing  the  measured  impedance  spectra  and  the 
model.  Since  the  impedance  of  the  DRT-model  does  not  contain  the 
differential  intercalation  capacity,  it  has  to  be  added  for  comparison 
with  the  measurement  data.  The  very  good  agreement  of  measured 
impedance  spectra  and  calculated  spectra  from  the  DRT  model  is 
shown  in  Fig.  8. 
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Fig.  8.  Measured  impedance  (points)  and  impedance  of  the  DRT  model  with  additional 
differential  intercalation  capacity  (straight  line).  Measurement  and  model  are  in  very 
good  agreement. 

4.2.  Evaluation  of  model  quality  in  time  domain 

The  simulation  result  of  the  validation  profile  with  a  number  of 
relaxation  times  of  N  =  100  is  depicted  in  Fig.  9.  A  good  agreement 
of  simulated  and  measured  battery  potential  is  achieved  covering 
the  complete  SOC  range.  For  validation  profiles  with  currents  up  to 
/max  =  1C  and  SOCs  down  to  5%,  the  potential  deviation  is  always 
lower  than  2%.  Even  though  there  is  a  visible  drift  of  the  simulated 
potential  over  time,  short-time  dynamics  are  modeled  with  high 
accuracy. 

For  the  profile  with  a  maximum  current  of  2C  the  overpotential 
of  the  model  is  consistently  too  high.  This  can  be  explained 
considering  the  measured  surface  temperature  of  the  cell.  Since  for 
the  validation  profile  of  2C  the  surface  temperature  increases  about 
5  K  (compare  Fig.  10),  the  true  impedance  of  the  system  is  no  longer 
represented  by  the  impedance  measured  under  isothermal  condi¬ 
tions.  Further,  the  thermal  conductivity  of  the  pouch  cell  perpen¬ 
dicular  to  the  electrodes  is  low.  Therefore  the  core  temperature  of 
the  cell  is  assumed  to  be  even  higher  than  the  measured  surface 
temperature.  This  higher  temperature  leads  to  a  decrease  of  the  cell 
impedance  and  hence  to  a  decrease  of  the  resulting  overpotential. 
In  order  to  account  for  this  effect  a  thermal  model  should  be 
coupled  to  predict  the  battery  temperature  and  adapt  the  modeled 
impedance. 

4.3.  Scalability 

The  aforementioned  scalability  is  given  by  the  freely  selectable 
number  of  relaxation  time  constants  N.  By  decreasing  N,  model 
complexity  and  computational  effort  is  decreased.  To  show  the 
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Fig.  7.  Distribution  of  relaxation  times  for  a  frequency  range  from  3.8  kHz  to  17  pHz.  The  SOC  is  varied  from  95%  to  5%  in  5%  decrements.  A  highly  SOC-dependent  process  can  be 
observed  between  1  Hz  and  10  Hz.  For  frequencies  higher  than  10  Hz,  an  SOC  dependence  is  not  observed,  as  can  be  seen  from  the  inset. 
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Fig.  9.  Validation  measurements  for  current  profiles  with  /max  ranging  from  C/4  to  2C.  During  the  validation  profile  with  /max,  the  cell  is  discharged  with  an  average  current  of  /max/2. 
The  straight  lines  show  the  simulated  cell  potential,  the  dotted  lines  the  measured  cell  potential. 
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Fig.  10.  Measured  surface  temperature  during  the  validation  profiles.  For  the  validation  profiles  with  /max  =  1C  and  2C  maximum  current,  a  significant  increase  of  the  surface 
temperature  can  be  observed. 
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Fig.  11.  Measured  impedance  at  an  SOC  of  50%  and  impedance  of  the  DRT-models  with 
the  orders  100,  20, 10  and  4. 


influence  of  the  number  of  relaxation  times  on  the  resulting  model 
impedance,  N  is  varied  with  N  =  [100  20  10  4].  The  impedance 
spectra  of  the  models  are  depicted  in  Fig.  11.  Naturally,  with 
decreasing  model  order  the  deviation  between  measured  and 
modeled  impedance  increases.  For  the  lowest  model  order  N  =  4 
the  semi  circles  of  the  RC  elements  are  clearly  distinguishable  and 
show  a  large  deviation  from  the  measured  impedance.  Flowever, 
already  for  N  =  20  the  model  coincides  with  the  measured 
impedance  very  well.  Besides  the  validation  in  frequency-domain, 
also  validation  in  time-domain  is  crucial.  As  Fig.  12a)  shows,  on 
a  scale  of  hours  the  model  quality  in  time  domain  does  not  differ 
significantly.  Flowever,  for  a  timescale  of  seconds  significant 
differences  are  observed  between  the  models  of  order  four  and  ten 
(compare  Fig.  12b)).  The  model  of  order  four  fails  to  reproduce  the 
measured  battery  voltage.  Among  the  remaining  models  of  higher 
order,  differences  are  smaller.  Flence  it  can  be  concluded,  that 
increasing  the  model  order  from  N  =  20  to  N  =  100  does  not  result 
in  a  significantly  higher  accuracy  while  the  computational  effort 
increases  with  a  factor  of  five.  In  sum,  the  model  of  order  ten  is  the 
best  tradeoff  between  accuracy  and  computational  effort  and  is 
therefore  recommended  for  the  modeled  frequency  range. 
Although  the  absolute  number  of  RC  elements  seems  to  be  rather 
high,  one  has  to  be  aware,  that  compared  to  other  approaches  in 
literature  [22],  the  covered  frequency  range  is  significantly  higher. 
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Fig.  12.  a)  Measured  and  simulated  battery  potential  for  a  validation  profile  with  /max  =  C/2.  Model  order  is  varied  from  100  down  to  4.  b)  Potential  difference  between  measured 
and  simulated  battery  potential  for  only  one  cycle.  The  model  of  order  four  fails  to  describe  the  dynamic  cell  behavior.  A  model  order  of  10  is  found  to  be  sufficient. 


Table  1 

Fraction  of  real  time  and  simulated  time  for  models  of  different  orders. 


Model  order,  N 

100 

20 

10 

4 

RC  elements/decade 

11.6 

2.3 

1.1 

0.5 

Wheal 

9% 

0.8% 

0.5% 

0.3% 

Considering  the  relative  number  of  RC  elements  per  frequency 
decade,  with  one  RC  element  per  decade,  this  value  is  rather  low. 

An  even  higher  model  quality  in  frequency  domain  could  be 
achieved  by  increasing  the  number  of  free  parameters  by  not 
presuming  a  certain  distribution  of  time  constants.  This  very 
interesting  option  was  implemented,  but  due  to  the  increase  of 
computational  complexity  a  full  automation  of  the  parameteriza¬ 
tion  processes  could  not  be  achieved  so  far.  Finding  a  solution  of 
this  problem  will  be  part  of  future  work. 

4.4.  Real-time  capability 

The  model  order  does  not  only  affect  accuracy,  it  also  affects  the 
computational  effort.  Table  1  shows  the  computation  effort  as  ratio 
of  real  time  and  simulated  time  Treai/’fsim  together  with  the  number 
of  RC  elements  per  frequency  decade  for  the  simulated  model 
orders.  All  simulations  were  performed  on  a  desktop  computer 
with  an  Intel®  Core™  i7  CPU  using  Matlab  R2009a  (Mathworks)  and 
PLECS  2.1  (plexim).  Except  for  the  model  of  order  100  all  simula¬ 
tions  could  be  performed  with  good  real  time  capability.  With 
respect  to  the  results  concerning  scalability,  the  model  of  order  10 
gives  the  best  trade-off  for  accuracy  and  computational  effort. 

5.  Conclusion 

The  introduced  model  is  parameterized  using  impedance  data 
from  a  Li-ion  pouch  cell  with  a  nominal  capacity  of  2  Ah.  It  has  the 
potential  for  application  in  on-board  diagnosis.  Because  of  the 
generalized  approach,  the  model  is  independent  of  material 
chemistry,  microstructure  and  capacity.  Therefore  no  a  priori 
knowledge  of  the  number  and  physical  origin  of  the  occurring 
processes  is  required.  Furthermore  it  combines  the  advantages  of 
automatable  parameterization,  scalability  and  a  structure  that  is 
directly  implementable  on  microcontrollers.  The  distribution  of 
relaxation  times  (DRT),  enables  (i)  a  direct  derivation  of  the  model 
from  measured  impedance  data  and  (ii)  a  full  automation  of  the 
parameterization  process.  Model  order  and  therefore  the 


computational  effort  are  scalable  via  the  discrete  number  of  time 
constants  considered  in  the  model.  Also,  the  considered  frequency 
range  can  be  easily  adjusted.  Therefore  the  model  is  versatile  and 
easily  adaptable  to  various  constraints. 

As  a  proof  of  this  concept,  the  linear  behavior  of  a  Li-ion  pouch 
cell  was  modeled  in  a  wide  frequency  range.  The  accuracy  in  time 
domain  was  validated  with  current  profiles,  which  show  that 
model  and  measurement  are  in  good  agreement  up  to  maximum 
currents  of  1C.  Furthermore  it  was  observed  that  for  the  investi¬ 
gated  high  power  cell,  even  currents  as  low  as  1C  cause  a  significant 
increase  of  the  battery  temperature.  As  this  temperature  influence 
is  assumed  to  be  responsible  for  most  of  the  prediction  error, 
coupling  of  an  additional  thermal  model  is  expected  to  be  highly 
beneficial. 
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